Non-Markovian effect on the geometric phase of a dissipative qubit 
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We study the geometric phase of a two-level atom coupled to an environment with Lorentzian 
spectral density. The non-Markovian effect on the geometric phase is explored analytically and 
numerically. In the weak coupling limit the lowest-order correction to the geometric phase is derived 
analytically and the general case is calculated numerically. It is found that the correction to the 
geometric phase is significantly large if the spectral width is small and in this case the non-Markovian 
dynamics has a significant impact to the geometric phase. When the spectral width increases, the 
correction to the geometric phase becomes negligible, which shows the robustness of the geometric 
phase to the environmental white noises. The result is significant to the quantum information 
processing based on the geometric phase. 

PACS numbers: 03.65.Vf, 03.65.Yz, 42.50.Lc 
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I. INTRODUCTION 

The concept of geometric phase (GP) in quantum sys- 
tem was originally introduced by Berry [l[ when he stud- 
ied the dynamics of a closed quantum system which un- 
dergoes an adiabatic cyclic evolution. He found that be- 
sides the usual dynamical phase, the system also acquires 
an additional phase which only depends on the geometry 
of the path traversed by the system during its adiabatic 
evolution. Since then this important notion has attracted 
much attention Q . It has been generalized in various as- 
pects, e.g. the GP for non-adiabatic evolution Q and 
for noncyclic evolution W. The GP has been observed 
experimentally in optical [5|, NMR 
ducting electronic circuit experiments 

Recently, the renewed interest in the investigation of 
GP comes from the application of GP to implement the 
logic gates in quantum computation . The purely ge- 
ometric nature of the phase makes such computation in- 
trinsically fault-tolerant and robust against certain types 
of classical fluctuation noise However, any realis- 

tic quantum system is inevitably coupled to its surround- 
ing environment, which would result in the loss of quan- 
tum coherence (i.e. the decoherence) of the quantum 
system itself and hence limit the implementation of the 
geometric quantum computation. Therefore, the study 
of the GP in open quantum systems becomes an impor- 
tant issue. For the GP of mixed states in open systems, 
Uhlmann was the first to make the attempt to define the 
mixed-state GP via state purification [13]. Sjoqvist et 
al. proposed an alternative definition for the nondegen- 
erate mixed-state density matrix under unitary evolution 
based on the interferometry (isj . This definition was fur- 
ther generalized to degenerate mixed state by Singh et 
al. [l9| and to the nonunitary evolution by Tong et al. 
[20| using the kinematic approach. The GP associated 



with such a nonunitary evolution is defined as (2C 
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where ek{t) and \ipk), respectively, are the eigenvalues 
and the eigenstates of the reduced density matrix of the 
quantum system, and T is the time after the system com- 
pletes a cyclic evolution when it is isolated from the envi- 
ronment. Taking the environment into account, the sys- 
tem no longer undergoes a cyclic evolution. Here a quasi- 
cyclic process with T — 27r/a;o: where wq is the frequency 
of the system, is considered in Eq. ([1]). Wang et al. 
defined a mixed-state GP in the context of Pantcharat- 
nam formula [2l] [ via mapping the density matrix to a 
nonunit vector ray in complex projective Hilbert space 
[2^ . The mixed-state GP has been observed in NMR 
system [2l,|2i|. 

Because the environment induced decoherence would 
affect the performance of the quantum computation using 
GP, the study of environment effects on the GP is highly 
desired. Many works along this line have been performed 
within Markovian approximation [25l - l30j . which is valid 
only when the interaction between the system and the 
environment is very weak and the environmental corre- 
lation time is very small. However, in many quantum 
information experiments, these conditions are not com- 
pletely satisfied. For example, in cavity QED experiment 
the imperfection of the cavity mirrors makes the cavity 
field having a Lorentzian spectrum expansion, which acts 
as an environment, would exert a strong non-Markovian 
effect on the atom in it 31]. There are also some works 
on the non-Markovian effect on the GP in dephasing en- 
vironments [33 - l36| . where there is no energy/information 
exchange between the system and its environment. 

In this work, we extend the study of the GP of open 
two-level system to the situation where there has an en- 
ergy/information exchange between the system and its 
environment. The environment is at zero-temperature 
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and has a Lorentzian spectral density, which corresponds 
to the radiation field as an environment being confined in 
a leaky cavity. We mainly concentrate on how the non- 
Markvoain effect affects the GP in different parameter 
regimes of the spectral density, and how the GP mani- 
fests its robustness against the decoherence. 

This paper is organized as follows. In Sec. |lTl we in- 
troduce the model of a qubit interacting with a vacuum 
environment and analyze its decoherence behavior. In 
Sec. mil we evaluate the GP of the qubit via performing 
analytical and numerical calculations. The correction ef- 
fect exerted by the environment on the GP in different 
parameter regimes of the Lorentzian spectral density is 
analyzed. Finally, a brief discussion and summary are 
given in Sec. IIVI 



II. THE MODEL 

Let us consider a system consisted of a two-level atom 
(qubit) coupled to a radiation field at zero temperature 
as an environment. The Hamiltonian of the system is 



H 



uJkalak + ^{gk(J+ak + h.c), (2) 



where a± and are the inversion operators and transi- 
tion frequency of the qubit, a^, and are the creation 
and annihilation operators of the k-th mode with fre- 
quency uJk of the radiation field, and gk is the coupling 
strength between the qubit and reservoir. Throughout 
this paper we assume h = 1. The system models the 
decoherence process of the atom via the amplitude de- 
caying under the Born-Markovian approximation, which 
results in the spontaneous emission of the two-level atom 
[3^ in quantum optics. The model is exactly solvable. 
The decoherence dynamics of the qubit is governed by 
the master equation [s^ 

p{t) = -zA(0[a+a_,p(t)]+r(i)[2a_p(i)a+ 

-a+a^p{t) - p{t)a+a-], (3) 



where the time-dependent parameters are given by 

AW = m . -Re||ll, (4) 

It is shown that c{t) satisfies 

c{t) + iujoc{t) + [ f{t~ T)c{T)dT = 0, (5) 



where f{t — t) ^ J J {uj)e~'^'^^*~'^^ duj is the environmen- 
tal correlation function with the spectral density defined 
as J (to) — X]fc ~ '^fc) t-^i^ initial condition 

c(0) = 1. The time-dependent parameters A(t) and T{t) 
play the roles of Lamb shifted frequency and decay rate of 
the qubit, respectively. The integro-differential equation 



([5]) contains the memory effect of the reservoir registered 
in the time-nonlocal kernel function and thus the dynam- 
ics of qubit displays non-Markovian effect. If the time- 
nonlocal kernel function is replaced by a time-local one, 
then Eq. ([3]) recovers the conventional master equation 
under Markovian approximation. 

It is obvious that the memory effect registered in the 
kernel function f{t — t) is essentially determined by the 
spectral density J{uj). In this work we explicitly consider 
that the spectral density has a Lorentzian form [soj 
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where W is the coupling constant between the qubit and 
the environment, and A defines the spectral width of the 
coupling at the resonance point wq. The Lorentzian spec- 
tral density describes that the vacuum radiation field as 
the environment is confined in a leaky cavity. Due to the 
leakage of the cavity field induced by the imperfection 
of the cavity mirrors, the spectrum of the cavity field 
displays a broadening at the resonance point regarding 
the atomic transition frequency ujq. In this case one can 
verify that the correlation function decays exponentially 
f{t~T) = iy2g-A(t-r) i3g^ which means that the pa- 
rameter A characterizes the correlation time of the en- 
vironment as Tc = A^^. If Tc is comparable with the 
typical time scale of the system, i.e. tq = l/wo, then 
the memory effect of the environment should not be ne- 
glected and the decoherence dynamics in this situation 
is non-Markovian. While ^ tq, the memory effect 
of the environment is negligible and the decoherence dy- 
namics is Markovian. In the ideal cavity limit A — > 0, we 
have lim^^o '-'^('^) = W"^ 5 {lo ~ ujt^) , which corresponds 
to a constant kernel f {t ~ t) — W"^ . Then the system 
reduces to the Jaynes-Cummings model with a vacuum 
Rabi frequency g = W. 

Going back to the general case, one can obtain the 
analytical form of c{t) by substituting the exponentially 
decaying f[t — T) into Eq. © as 



c(t) 



(A + i^o)t 



cosh(y) + -smh(y) 



(7) 



where = VA^ — 4W'^. So the parameters in the master 
equation can be calculated readily 



2W^ 



(8) 



One can see from Eqs. (|S]) when A is much larger than 
other frequency scale, the decay rate tends to a constant 
value Fq = 2W'^/X, which just characterizes the decoher- 
ence behavior of the qubit under the Markovian dynam- 
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III. GEOMETRIC PHASE CORRECTED BY 
THE ENVIRONMENT 

A. Analytical analysis 

In the following we compute explicitly the GP of the 
qubit. We assume that the initial state of the qubit is 
chosen as 




\m) 
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(9) 



where |+) and |— ) are the excited and ground states 
of the qubit, respectively. This state corresponds to a 
vector in Bloch sphere with polar angle ^o- The time- 
dependent reduced density matrix of the qubit under the 
initial condition ([9|) can be obtained straightforwardly 
from the master equation ([3]) 



Pit) 



COs2|i|c(i)|2 
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(10) 



It is noted that besides the off-diagonal elements, the di- 
agonal elements of the reduced density matrix p{t) also 
change with time in our model. It is just this time- 
dependence of the diagonal elements of p{t) character- 
izing the energy exchange between the qubit and its en- 
vironment that makes our system shows dramatic differ- 
ence to the dephasing model i32i-t36i] . To calculate the 
GP of the qubit, we must firstly get the eigensolution of 
the reduced density matrix ((TU]). The eigenvalues of the 
above reduced density are readily calculated. 



e±it) 



1± W|c(i)|2sin'6lo-l- 2|c(t)|2cos 



1 



(11) 

It is obvious that the eigenvalue £-(0) = 0, which, from 
Eq. ([TJ, means that the component of the eigenstate 
corresponding to the eigenvalue £_ gives no contribution 
to the GP. Thus we only need to consider the eigenstate 
corresponding to the eigenvalue e+ 



-iuQt 



COS Q I 



sin Q I 



(12) 



where 



COS 6 = 



2{\c{t)\ 



2 C0S2 ^ 



,{t)\^ 



sm Oq 



|c(t)|2cos2f 



(13) 

Below we calculate the GP. Eq. ([8]) shows that the 
frequency shift of the qubit is zero for the Lorentzian 
spectral density [38l |. So the period of the environment 
disturbed atom is the same as the one for a bare atom. 
Then the GP of the qubit acquired after a period T = 
27r/a;o can be calculated as 



uiQ cos 



(14) 



FIG. 1: (Color online) Numerical result of the GP as a func- 
tion of the coupling constant and the initial angle, (a) Marko- 
vian dynamics with A = S.Oojq and (b) non-Markovian dynam- 
ics with A = O.OSwo. 



Here we point out that the expression of the GP (fT4|) 
has been obtained by using the kinematic approach 20j . 
However, it can be shown that the expression is ex- 
actly the same as the one obtained by the definition 
of Pantcharatnam's phase of a nonunit vector ray [2^ . 
which exhibits some essential features of the GP in such 
systems. 

From Eq. ((HI), one can notes that the effect of the 
environment enters into the GP solely via the time- 
dependent factor |c(t)|2 of the excited-state population. 
When the environment is absent, then |c(t)|2 = 1 and 

the GP reduces to = tt (1 -f cos^o), which is just 
the GP acquired by a two-level atom in the unitary dy- 
namics. Depending on the decoherence dynamics being 
Markovian or non-Markovian, the short-time dynamics 
of |c(t)|2 shows remarkably different behaviors. In the 
Markovian dynamics, |c(i)|2 decays monotonically and 
finally approaches zero. In this case, the larger the decay 
rate Pq is, the larger the correction of the GP should be. 
In the non-Markovian dynamics, contributed from the 
memory effect of the environment |c(t)|2 shows transient 
oscillation with time, which naturally induces a correc- 
tion of the GP different to the one in Markovian dynam- 
ics, as shown in the following. 

Up to the first order of coupling strength 11^2^ j the 
weak coupling limit, we have 



$(1) ^ $(0) 

a 9 



■ 2. 

— sm t^o 



cos Uq 
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(15) 

Tra;)] . Besides 



where z (x) = 

the leading term $1,'^'' corresponding to the well-known 
GP acquired under the unitary dynamics, the second 
term, which is quadratic in W/ujq with a A-dependent 
coefficient z{\/ujq), is the lowest-order correction to the 
GP induced by non-unitary dynamics due to the inter- 
action with the environment. It is easy to check that 
z{X/u!q) is a monotonically decreasing function with the 
increase of A, so we can expect that the GP shows a 
larger deviation to for a small A than a large one. In 
particular, in the ideal cavity limit A — > 0, the function 
z{X/ujq) arrives at its maximum Att^ /3, where the GP has 
a largest correction in this weak coupling (or small W) 
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regime. On the other hand, when A ^ wq, W, one can 
verify z{X/ujq) — >■ 2Tr'^ujQ/X. Consequently, the GP in this 
Markovian hmit is 

<I'^^)=<i>r-^-^^o(l + ^^), (16) 

which shows very similar form to the result of (26| ob- 
tained from the Born-Markovian master equation. 
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FIG. 2: (Color online) The exact GP and $g as a function of 
the coupling constant under different spectral width A. The 
initial polar angle is taken as 9o = n/3. 
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FIG. 3: (Color online) The exact GP and as a function of 
the initial polar angle under different spectral width A. The 
coupling constant is = 0.2a;o. 



B. Numerical results 

In the following we give the exact GP via numerically 
evaluating the integration in Eq. (|14l) . In Fig. [T] we 
plot the exact GP as a function of the coupling constant 



and the initial polar angle. When A = S.Owqi the cor- 
relation time of the environment Tc — 0.2/ujo is much 
less than the typical time scale of the system tq. In 
this situation, the memory effect of the environment is 
negligible and the dynamics of the system is Markovian. 
Consequently, the energy and the information flow single- 
directionally from the system to the environment, which 
results in the dissipation to the system. The dissipation 
time scale in this Markovian dynamics is characterized 
by TD = l/To = X/2W^. We can see in Fig. [BJa) that 
the GP in this situation decreases monotonically with the 
increase of the coupling constant. This is understandable 
that a larger coupling constant would induce a stronger 
dissipation and a short decoherence time scale td to the 
qubit system. With the increase of W, td approaches 
T, which means that the dissipation becomes more and 
more notable within the time scale T. On the other hand, 
when W is small the GP shows small deviation to the 
unitary one <I>g°' . In this situation, td is larger than T , 
making the dissipation negligible within the time scale 
T during which the GP is accumulated. This interesting 
phenomenon just reflects the resilient ability of the GP 
to the environment's disturbance, especially in weak cou- 
pling limit. A similar phenomenon is also observed in the 
dephasing model [s^. When A = O.OScjq, Tc ^ tq and 
the non-Markovian effect induced by the memory effect 
of the environment to the system would show its dis- 
tinct impact on the the dynamics of the system. Besides 
the dissipation, the environment also exerts a dynami- 
cal backaction on the system The backaction effect 
is reflected in that the energy and the information flow 
back and forth between the system and the environment. 
In this situation we can see in Fig. [IJb) that the GP 
oscillates with the increase of W and then tends to a def- 
inite value, which is qualitatively different to the above 
Markovian situation. The oscillation is just the manifes- 
tation of the dynamical backaction. It is very clear that 
the non-zero character of the GP in the case of large 
W is due to the counteracted competition between the 
backaction and dissipation exerted by the non-Markovian 
environment partially, which weakens the decoherence of 
the qubit system. 

A more clear comparison of the GP with different A as 
a function of W when 9o = tt/S is shown in Fig. [5] As 
comparisons, the perturbed results for A = and S.Owq 
are also presented in this figure and they agree well with 
the corresponding exact one only in the weak coupling 
limit. We can see obviously that in the weak coupling 
regime the environment with a large A induces a smaller 
correction to the GP than the one with a small A. That is, 
the non-Markovian effect has a strong correction to the 
GP in this weak coupling regime. This is consistent with 
the result obtained in Ref. [4l|, where a phenomeno- 
logical analysis of the non-Markovian effect on the GP 
is performed. While in the strong coupling regime the 
situation is opposite, because the system with a large-A 
shows stronger decoherence than the one with a small A 
where the non-Markovian effect has a strong correction 
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to the GP. 

Even till today it is still difficult to access experimen- 
tally the strong coupling regime in the cavity QED plat- 
form. For example, in a recent microtoroidal resonator 
experiment J2j , the achieved coupling strength between 
the atom and the cavity field is 90MHz, while the damp- 
ing rate of the cavity field is about 180MHz. Such a 
large damping rate would induce a noticeable spectrum 
expansion of the cavity mode. It is of course interesting 
to examine the GP in this bad-cavity and weak coupling 
regime. In Fig. [3] we plot the GP as a function of the 
initial polar angle in the weak coupling regime for dif- 
ferent spectrum width A. One can see that the GP ap- 
proach more and more closely the unitary GP $g'^'* with 
the increase of A. This accounts for once again that the 
GP shows more strong resilient ability to the noise of 
the Markovian environment with a large A. Since this 
resilient character is obtained in weak coupling and bad- 
cavity condition, we argue that it is accessible by modern 
cavity QED experiment (i^ . 

IV. SUMMARY AND DISCUSSION 

We have studied exactly the GP for a qubit in an 
amplitude decaying environment using the kinematic ap- 
proach. We have evaluated the non-Markovian effect on 
the GP. It is demonstrated that the non-Markovian en- 
vironment has a more significant correction to the GP in 
the weak coupling regime due to the strong short-time 
correlation between the qubit and the environment. In- 
terestingly, our result also indicates the insensitivity of 
the GP to the Lorentzian environment in the Markovian 
regime when the coupling is weak. It just manifests the 
resilient ability of the GP to the environmental noise. 



This results elucidate that the GP in this cavity QED sys- 
tem is fault-tolerant not only against the classical noise 
induced by the parameter fluctuation [ll|, [11], but also 
against the quantum noise. This has significant meaning 
in geometric quantum information processing. 

Our model, as a basic model of quantum optics, is 
particularly relevant to the cavity QED experiments. In 
this respect some remarkable experiments, such as the 
efficient coupling of the trapped atoms with cavity field 
[43I [4^ . have been performed successfully. Practically, 
any phase variation is observed only via some kind of in- 
terferometry between the involved state and certain se- 
lected reference state. For example, the GP for mixed 
state has been observed via designing a quantum network 
in NMR system to realize the interferometry [23I, [2^ . 
This provides a clue to observe GP in cavity QED sys- 
tem. Although the GP has not been observed in the 
cavity QED system, a quantum network using the recent 
developed microtoroidal resonator [42] has been proposed 
jisj based on the input-output process of photons p6j . 
If an effective interferometry could be realized in this 
quantum network, then the GP would be expected to 
be observable in cavity QED system. Our work on the 
assessment of the environmental effect on the GP, espe- 
cially in non-Markovian regime, is of great importance 
in using the GP in cavity QED system to implement the 
quantum gates. 
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